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Abstract 

We compute the Hamiltonian for spherically symmetric scalar field collapse 
in Einstein-Gauss-Bonnet gravity in D dimensions using slicings that are regular 
across future horizons. We first reduce the Lagrangian to two dimensions using 
spherical symmetry. We then show that choosing the spatial coordinate to be a 
function of the areal radius leads to a relatively simple Hamiltonian constraint 
whose gravitational part is the gradient of the generalized mass function. Next 
we complete the gauge fixing such that the metric is the Einstein-Gauss-Bonnet 
generalization of non-static Painleve-Gullstrand coordinates. Finally, we derive 
the resultant reduced equations of motion for the scalar field. These equations 
are suitable for use in numerical simulations of spherically symmetric scalar field 
collapse in Gauss-Bonnet gravity and can readily be generalized to other matter 
fields minimally coupled to gravity. 
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1 Introduction 

The rising popularity of higher dimensional gravitational theories such as string theory 
calls for the consideration of higher order curvature terms in the action [Tj. If we require 
the action to yield covariant equations that are at most second order in derivatives of the 
metric, then we must modify the action by adding Lovelock polynomials [2], as these are 
the only known combinations of higher order curvature terms that satisfy this criterion. 
These terms are identically zero in four dimensions but are relevant for D > 4. The 
second order Lovelock polynomial is known as the Gauss-Bonnet (GB) term. 

A considerable amount is known about static black hole solutions in Lovelock gravity 
[3] and there have been several analytic studies of the dynamics of black hole horizons in 
Einstein- Gauss-Bonnet (EGB) gravity. For example, the effect of GB terms on the end 
state of gravitational collapse has been investigated in generalized Vaidya co-ordinates 
As well, dynamical horizons were studied in the context of EGB gravity in [6J and 
Surprisingly, comparatively little is known about the formation of black holes in 
EGB gravity. This situation stands in strong contrast to that in general relativity, where 
numerical studies have yielded much information about the formation of black holes from 
gravitational collapse Here we take the first steps in this direction by deriving the 
Hamiltonian equations describing gravitational collapse of a spherically symmetric scalar 
field in Einstein- Gauss-Bonnet (EGB) gravity. 

The general form of the Hamiltonian for Lovelock gravity was first derived in [9] using 
a covariant formalism. This work did not write down the explicit form of the Hamilto- 
nian in terms of phase space variables and was therefore not of much use practically 
for investigating collapse. Further progress came when the Hamiltonian formulation of 
vacuum EGB in D = 5 in terms of Kuchar's geometrodynamic formalism was employed 
to investigate black hole thermodynamics p^. More recently, the role of surface terms 
and conserved charges in EGB was explored fTIj [T3] and the mass-energy of a EGB space 
time was constructed [H]. 

In the present work, we present a complete Hamiltonian analysis of a massless scalar 
field minimally coupled to EGB gravity. We first derive a partially reduced Hamiltonian 
in a particular family of gauges in which a function of the area radius is used as the spatial 
coordinate. This yield slicings that are regular across future (or past) horizons. Moreover, 
the Hamiltonian constraint can be simply expressed in terms of the generalized mass 
function first derived in [H]. We then fully reduce the theory in generalized Painleve- 
GuUstrand (PG) coordinates, obtaining the reduced equations of motion for the scalar 
field and its conjugate. The reduced equations are in a form that is useful for the 
analysis of gravitational collapse up to and beyond horizon formation either analytically 
or numerically as done for Einstein gravity in |8j. Our analysis can readily adapted to 
other types of minimally coupled matter fields. 

The paper is organized as follows. In Section [2] we dimensionally reduce the spherically 
symmetric, D- dimensional action to one space and one time co-ordinate. We then write 
the action in a form that manifestly gives second order derivatives of the metric in the 
covariant equations of motion. Section 12] presents the Hamiltonian analysis and derives 
the partially gauge fixed equations. Finally, in Section H] we complete the gauge fixing 
and calculate Hamilton's equations of motion for a massless scalar field in the presence 



2 The Action 



of gravity in generalized PG co-ordinates. We then conclude in section O Some technical 
details are relegated to an Appendix. 



2 The Action 

We start with the Einstein-Hilbert action with a Gauss Bonnet term minimally coupled 
to a massless scalar field, 

S = Sg + Sm = Seh + Sgb + Sm-, (2.1) 

where Sq is the gravitational part the action, with contributions from the Einstein-Hilbert 
term, Seh, and the Gauss-Bonnet term, Sgb and Sm comes from the scalar field. These 
terms are given by 
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and 

Sm = -\j d''x^-~g\y^\\ (2.4) 

where Gr, is Newton's gravitational constant, D is the number of dimensions, g is the 
determinant of the metric, R is the Ricci scalar, R^^ is the Ricci tensor, R^ypa is the 
Riemann tensor, a is the Gauss Bonnet coupling constant and is the scalar field. 
Inserting the spherically symmetric ansatz 

ds^ = gab{t, x)dx''dx^ + r\jdx'dx^ (2.5) 

where 'jij is the metric of the unit n sphere, x and t are spatial and time co-ordinates, r 
is the areal radius and the indices a and b now run from to 1. Performing the angular 
integrals in the action, we get 
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3 Hamiltonian Analysis 



where n = D — 2 is the number of angular dimensions and the surface area of an 

n sphere, g and R are the determinant and Ricci scalar, respectively, of the 2D metric, 

9ab 

We now consider only the gravitational part of the action. Integrating by parts and 
defining = r"/^ ^16j , we find 
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A key relationship in two dimensions is the following: 

dxdty^i^ [-\V(I)\^R + 2(n0)2 - 2(VaV60)(V'^VV)] = 0, 
This allows for further simplification of the action 
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Equation (12. lip is now in a form that clearly yields Euler-Lagrange equations that are 
second order in derivatives of the metric, as expected from Gauss-Bonnet theory. 



3 Hamiltonian Analysis 

We work with the dimensionally reduced two-geometry and write Qab in ADM form 



ds'^ = {-a'^dt^ + {dx + Ndtf) 



(3.1) 



where p, a and are functions of x and t. Since the theory is diffeomorphism invariant 
we know that the total Hamiltonian will be of the form [S] 



H = j dx{ag + NJ^) 
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3.1 Gravitational Hamiltonian 



The metric functions a and are Lagrange multipliers and Q and J-" are the Hamiltonian 
and diffeomorphism constraints. Moreover, as noted in (TU] the form of J-" is completely 
specified by the choice of parametrization and the fact that J-" generates spatial diffeo- 
morphisms. Specifically, in the present case the gravitational part is: 

j^G = Hpp' - n; + n^(p' (3.3) 

We first focus on the gravitational sector of the theory. 



3.1 Gravitational Hamiltonian 

The Hamiltonian analysis for vacuum 5-D Gauss-Bonnet theory was done in [10] using 
the geometrodynamical formulation of Kuchar [11]. The algebra gets considerably more 
complicated in higher dimensions and as additional orders in derivatives are added in 
the more general case of Lovelock gravity. Here we take a simpler approach that is well 
adapted to the class of gauges we wish to consider. 

First of all it is important to note that for the action Sg only the term containing 
the Ricci scalar contains derivatives of the metric. Moreover, since the 2-D Ricci scalar 
for the parametrization (13.11) is 

R=-2ap (3.4) 

where □ refers to the D'Alembertian of the metric: {^—a'^dt^ + {dx + Ndt)'^), it follows 
that after an integration by parts, the action contains terms at most linear in the time 
derivative of p. This in turn implies that the gravitational part of the action can formally 
be written in the form: 

Sg = j dtdx(^BS,^,p,a,N)p + Bo{^,^,p,a,N)'^ . (3.5) 

Thus 

H, = 5i(0,0,p,a,iV) (3.6) 
which in principle at least can be solved for 0: 

= 0(Hp,p,0,(T,iV) (3.7) 

Note that in the above and what follows we do not specify separately dependence on 
spatial derivatives of a field separately from dependence on the field itself. 
The Legendre transformation then yields a Hamiltonian of the form: 

Hg = j dx{U^^{Up, p, 0, a, N) - 5o(0(H„ p, 0, a, N),p, 0, a, N)). (3.8) 

where as indicated is an implicit function of the other phase space variables as obtained 
from (13. 6 p but does not depend on H,^. The crucial property of the above is that the p 
terms have cancelled and the only dependence on H,^ is in the first term. 

We now make the above explicit for the action Sg = Seh + Sgb- In particular, we 
have: 

Sg = Q J dxdtiA^^-^ + ^30=^ + ^20^ + Ai0 + Ao) (3.9) 



3.2 Partial Gauge Fixing 



where = 8(n - l)z/(")/(167rGi5n), 
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Here 7 is defined as 7 = P i6(n-i)(Ti.-2) ^j^^ ^^^^ primes represent differentiation 
witli respect to tfie t and x coordinates respectively. 

For completeness we note that the conjugate momentum of is: 

= Ai + 2A2(j) + 3^30' + 4 A4(/)^ (3.15) 



3.2 Partial Gauge Fixing 

The above general discussion suggests a particularly convenient family of gauges, namely 
those for which the areal radius is time independent. We therefore consider the class of 
gauge fixing conditions: 

X/ = 0-/(x)^O (3.16) 



3.2 Partial Gauge Fixing 



If f{x) = x^/", then this gives x = r, i.e. the spatial coordinate is equal to the areal 
radius, but for the moment we will consider the more general case. 

The ~ denotes the fact that the gauge fixing condition can only be imposed on initial 
data. If it is a good gauge fixing condition then the requirement that it be preserved by 
the time evolution will yield a condition on one of the lagrange multipliers a or N. Xf 
must be second class with some linear combination of the constraints. Clearly: 
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From Eq.f l3.6p we can deduce that the required condition on the phase space variables is: 
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which determines the shift N in terms of a, p, lip and cj). The latter is no longer dynamical. 

In order to set x strongly equal to zero in the Hamiltonian (as opposed to the equations 
of motion) we must use the formalism of Dirac brackets. Fortunately, since is no longer 
a phase space variable, the Poisson bracket of all phase space variables with x is zero. 
This means that the Dirac bracket is the same as the Poisson bracket. 

The next step in the gauge fixing procedure is to use the diffeomorphism constraint 
to solve for II,;!, in terms of the other phase space variables. While this is indeed possible, 
it is not necessary in this case because the general form of the Hamiltonian (13. 8 p implies 
that when = 0, II,^ is not present in the Hamiltonian. Thus in this family of gauges, 
the gravitational part of the Hamiltonian is simply: 



Hg = - I dxBn. 
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where it is understood that in Bq, the shift is replaced by its value as determined by 

(EUD. 

To find the Hamiltonian with = we use fl3.19p . noting that -Bo(0 = 0) = Aq{p = 0) . 
This gives 
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where Hp is calculated to be 
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3.2 Partial Gauge Fixing 



which is a cubic equation that can be solved for N/a. As noted in [TU] there can in 
principle be multiple real roots, depending on the parameters involved, which would 
yield different dynamical spacetimes. In that work, the boundary conditions effectively 
resolved this potential ambiguity. In the appendix we show that in our case as well there 
is a unique real branch, and hence unique dynamics, for the fiat slice, generalized P-G 
coordinates that we consider. 

We therefore assume that there is a unique solution that specifies implicitly N/ cr as a 
function of the remaining phase space variables and we use this to write the Hamiltonian 
in terms of one Lagrange multiplier, cr, and one constraint 
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where V is the coefficient of R in the action, equation (12.111) iV,^ = 7^0^ + 8(^rxy' 
and g^^ is given by 



After some algebra, the Hamiltonian can be written as 
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where we have defined the mass function Ai [1T| as: 
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Using: 
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one can verify that M. agrees with the mass function defined in equation (2.7) of [14J. 
It is therefore the EGB generalization of the Misner-Sharp mass function. It is invariant 
under coordinate transformations that preserve spherical symmetry and asymptotes to 
the ADM mass for asymptotically fiat spacetimes. Clearly Ai is constant for vacuum 
spacetimes. 



3.3 Adding the Matter Contribution 



Eq.f l3.25p is a quadratic that can be inverted to express (g^^) and (via fl3.23p ) N/a in 
terms of the mass function: 



l±Jl + 167(72 - l)20-2-2/n>f /nQ (3.27) 



From the above, we see that in vacuum, when = M is constant, the above yields the 
two branches of the well-known vacuum solution of EGB theory. For the following, we 
choose the plus-branch in f l3.27p in order that our solutions reduce to those in Einstein 
gravity in the limit 7 — )■ 0. However, nothing in principle precludes the other choice, 
which would have interesting consequences for the end-point of gravitational collapse. 
This in turn puts the ADM metric in the form: 

ds^ = e^" [-a^e^'{g^^)de ± 2a^l - e^P{g^^)dtdr + cir^) . (3.28) 

with {g^^) given in fl3.27p above. Thus, the partially gauge fixed geometry is now implic- 
itly specified in terms of p, lip and the the lapse a. 

Finally, noting that the condition for the existence of an apparent horizon is that the 
null expansion vanishes: 

|Vr|2 = (7ii(r')' = (3.29) 

it is easy to verify that by choosing the -f- (-) in front of the off-diagonal term fl3.28p one 
can guarantee that the slicings are regular across future (past) horizons. The exception to 
this is the gauge choice = (Schwarzschild gauge) which breaks down at all horizons. 

3.3 Adding the Matter Contribution 

In order to add the scalar field contribution to the action we dimensionally reduce equa- 
tion ([MD to get 

Sm = TT- / dxdt^/^(j)^\V^\^ = (3.30) 
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From this we can write the matter action as 
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From equation (13.321) it is clear that the total partially gauge-fixed Hamiltonian is 
given by 
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where (^) is understood to the be real branch of the solution to (13.211) . The Hamiltonian 
(I3.33P is the generalization to G-B gravity of the partially reduced Hamiltonian derived 
in [15]. Note that it is valid for any gauge fixing condition Xf the form (I3.16p . For 
concreteness we henceforth specify: 
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which chooses the areal radius as the spatial coordinate: x = r = </)^/". 

At this stage, the dynamical system consists of four phase space variables (p, Hp, xjj, H^) 
and there is a single first class constraint, namely the Hamiltonian constraint. We now 
look at the complete reduction to only the physical degrees of freedom associated with 
the scalar field. 



4 Completely Reduced Theory 

For our second gauge choice we follow [8] and choose 
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It is clear from (13. ip and (12. 5 p that this choices yields spatial slices that are fiat. It 
therefore generalizes P-G coordinates. To further understand the significance of this 
choice, note that using (I3.27P in (I3.28P and setting x to zero in the limit as 7 — > 
(choosing the plus sign), we get g^^ — j- 1 — 2GM./r. Inserting this into (I3.28P and setting 
^ to zero we find that the metric becomes 



ds' 



2 / 2GM\ , 2 o 
-aM 1 — j dt^ + 2a 



2GM 



■dtdr + dr"^. 



(4.2) 



This shows that the gauge choice ^ = puts the metric into the Gauss-Bonnet, non- 
static generalization of the Painleve Gullstrand co-ordinates. Using ^ = also allows us 
to write the consistency condition for the first gauge choice, equation (I3.27p . as 
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5 Conclusion 



To ensure that the second gauge condition is conserved in time we must insist that 
{d/dt){e^p - 1) = {{e^P -1),H} = 0^ iS/6Ilp)H = 0. From flX^ we know that the Up 
dependence in the Hamiltonian is in the {N/a) terms. So we can write 
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where it is understood that we use (14.31) to write {N/a) in terms of the mass function 
M. 

Using X = ^ = we can now write down the fully gauge fixed Hamiltonian as 
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Using Hamilton's equations we find 



(4.6) 



ip = (J 



+ 



H 



H 



(4.7) 



(4.8) 



These equations, along with the consistency conditions (14. 3 p and (14. 5 p and the Hamilto- 
nian constraint 
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determine the evolution of a collapsing scalar field. 
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5 Conclusion 

Lovelock gravity has become increasingly popularized in recent years, as other theories 
have required higher dimensional analogues to General Relativity. We have presented 
here the Hamiltonian formulation of the lowest order non-trivial Lovelock theory, Gauss- 
Bonnet gravity, minimally coupled to a scalar field in the spherically symmetric case. 

The reduced Hamiltonian system (14. 51 - 14. 9]) is well suited for the study of gravitational 
collapse of a scalar field along with the resultant critical phenomena. Our calculations 
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utilize a gauge choice that selects non-static Painleve-GuUstrand coordinates, a coordi- 
nate system that is regular across the horizon during black hole formation. Moreover, this 
system is easily generalized to any matter field whose action does not contain derivatives 
of the metric. 

Future work will involve numerically studying gravitational collapse of a scalar field 
in this scenario to see how non-linearity in curvature modifies the well-known critical 
phenomena [1] in Einstein gravity. 
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7 Appendix 

We now solve equation 13.211 for N/a. Equation 13.211 can be written in the form 

c,(^)\c,[^)-n,^o (7.i; 

where 



Providing C1/C3 > 0, there is only one real solution, namely: 

(lOSn, + 12V3Cl^4Cf/Cs + 27n2]'^' 2Ci 



^'^^ f lOSn^ + 12V3ClJACf/Cs + 27Ul' 
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The condition that C1/C3 > can be simplified to 

^ 7n2</)2("-2)/" + ^\ - 1 > (7.4) 



470'^ V ^ - 1 

for the solution to equation 17.11 to be real. This won't be identically true so to proceed 
we must be careful with our gauge choices or possibly our allowed range of 7. For the 
fiat slice coordinates of Section (jl]) x = r (ie 0' = (r;,/2)0^~^/") and e^^ = 1 so that the 
condition is in fact identically satisfied. 
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Abstract 

We compute the Hamiltonian for spherically symmetric scalar field collapse 
in Einstein-Gauss-Bonnet gravity in D dimensions using slicings that are regular 
across future horizons. We first reduce the Lagrangian to two dimensions using 
spherical symmetry. We then show that choosing the spatial coordinate to be a 
function of the areal radius leads to a relatively simple Hamiltonian constraint 
whose gravitational part is the gradient of the generalized mass function. Next 
we complete the gauge fixing such that the metric is the Einstein-Gauss-Bonnet 
generalization of non-static Painleve-Gullstrand coordinates. Finally, we derive 
the resultant reduced equations of motion for the scalar field. These equations 
are suitable for use in numerical simulations of spherically symmetric scalar field 
collapse in Gauss-Bonnet gravity and can readily be generalized to other matter 
fields minimally coupled to gravity. 



1 Introduction 



1 Introduction 

The rising popularity of higher dimensional gravitational theories such as string theory 
calls for the consideration of higher order curvature terms in the action [Tj. If we require 
the action to yield covariant equations that are at most second order in derivatives of the 
metric, then we must modify the action by adding Lovelock polynomials [2], as these are 
the only known combinations of higher order curvature terms that satisfy this criterion. 
These terms are identically zero in four dimensions but are relevant for D > 4. The 
second order Lovelock polynomial is known as the Gauss-Bonnet (GB) term. 

A considerable amount is known about static black hole solutions in Lovelock gravity 
[3] and there have been several analytic studies of the dynamics of black hole horizons in 
Einstein- Gauss-Bonnet (EGB) gravity. For example, the effect of GB terms on the end 
state of gravitational collapse has been investigated in generalized Vaidya co-ordinates 
As well, dynamical horizons were studied in the context of EGB gravity in [6J and 
Surprisingly, comparatively little is known about the formation of black holes in 
EGB gravity. This situation stands in strong contrast to that in general relativity, where 
numerical studies have yielded much information about the formation of black holes from 
gravitational collapse Here we take the first steps in this direction by deriving the 
Hamiltonian equations describing gravitational collapse of a spherically symmetric scalar 
field in Einstein- Gauss-Bonnet (EGB) gravity. 

The general form of the Hamiltonian for Lovelock gravity was first derived in [9] using 
a covariant formalism. This work did not write down the explicit form of the Hamilto- 
nian in terms of phase space variables and was therefore not of much use practically 
for investigating collapse. Further progress came when the Hamiltonian formulation of 
vacuum EGB in D = 5 in terms of Kuchar's geometrodynamic formalism was employed 
to investigate black hole thermodynamics p^. More recently, the role of surface terms 
and conserved charges in EGB was explored fTIj [T3] and the mass-energy of a EGB space 
time was constructed [H]. 

In the present work, we present a complete Hamiltonian analysis of a massless scalar 
field minimally coupled to EGB gravity. We first derive a partially reduced Hamiltonian 
in a particular family of gauges in which a function of the area radius is used as the spatial 
coordinate. This yield slicings that are regular across future (or past) horizons. Moreover, 
the Hamiltonian constraint can be simply expressed in terms of the generalized mass 
function first derived in [H]. We then fully reduce the theory in generalized Painleve- 
GuUstrand (PG) coordinates, obtaining the reduced equations of motion for the scalar 
field and its conjugate. The reduced equations are in a form that is useful for the 
analysis of gravitational collapse up to and beyond horizon formation either analytically 
or numerically as done for Einstein gravity in |8j. Our analysis can readily adapted to 
other types of minimally coupled matter fields. 

The paper is organized as follows. In Section [2] we dimensionally reduce the spherically 
symmetric, D- dimensional action to one space and one time co-ordinate. We then write 
the action in a form that manifestly gives second order derivatives of the metric in the 
covariant equations of motion. Section 12] presents the Hamiltonian analysis and derives 
the partially gauge fixed equations. Finally, in Section H] we complete the gauge fixing 
and calculate Hamilton's equations of motion for a massless scalar field in the presence 



2 The Action 



of gravity in generalized PG co-ordinates. We then conclude in section O Some technical 
details are relegated to an Appendix. 



2 The Action 

We start with the Einstein-Hilbert action with a Gauss Bonnet term minimally coupled 
to a massless scalar field, 

S = Sg + Sm = Seh + Sgb + Sm-, (2.1) 

where Sq is the gravitational part the action, with contributions from the Einstein-Hilbert 
term, Seh, and the Gauss-Bonnet term, Sgb and Sm comes from the scalar field. These 
terms are given by 
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d^xy/-gR, (2.2) 



SGB = aJ d^xv^ (i?" - AR^.R'^'' + R^.p^R'^"'"'^ , (2.3) 



and 

Sm = -\j d''x^-~g\y^\\ (2.4) 

where Gr, is Newton's gravitational constant, D is the number of dimensions, g is the 
determinant of the metric, R is the Ricci scalar, R^^ is the Ricci tensor, R^ypa is the 
Riemann tensor, a is the Gauss Bonnet coupling constant and is the scalar field. 
Inserting the spherically symmetric ansatz 

ds^ = gab{t, x)dx''dx^ + r\jdx'dx^ (2.5) 

where 'jij is the metric of the unit n sphere, x and t are spatial and time co-ordinates, r 
is the areal radius and the indices a and b now run from to 1. Performing the angular 
integrals in the action, we get 
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3 Hamiltonian Analysis 



where n = D — 2 is the number of angular dimensions and the surface area of an 

n sphere, g and R are the determinant and Ricci scalar, respectively, of the 2D metric, 

9ab 

We now consider only the gravitational part of the action. Integrating by parts and 
defining = r"/^ ^16j , we find 
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A key relationship in two dimensions is the following: 
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which is true up to a total divergence. This allows for further simplification of the action 
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3 Hamiltonian Analysis 

We work with the dimensionally reduced two-geometry and write gab in ADM form 

ds^ = {-a^de + {dx + Ndtf) , (3. 1) 

where p, a and are functions of x and t. Since the theory is diffeomorphism invariant 
we know that the total Hamiltonian will be of the form [U] 



H = dx{(rg + NJ^) 



(3.2) 



The metric functions cr and A^ are Lagrange multipliers and Q and J-" are the Hamiltonian 
and diffeomorphism constraints. Moreover, as noted in the form of J-" is completely 
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specified by the choice of parametrization and the fact that T generates spatial diffeo- 
morphisms. Specifically, in the present case the gravitational part is: 

Tg = Hpp' - n; + n<^0' (3.3) 

We first focus on the gravitational sector of the theory. 
3.1 Gravitational Hamiltonian 

The Hamiltonian analysis for vacuum 5-D Gauss-Bonnet theory was done in [10] using 
the geometrodynamical formulation of Kuchar [TT]. The algebra gets considerably more 
complicated in higher dimensions and as additional orders in derivatives are added in 
the more general case of Lovelock gravity. Here we take a simpler approach that is well 
adapted to the class of gauges we wish to consider. 

First of all it is important to note that for the action 5*^ only the term containing 
the Ricci scalar contains derivatives of the metric. Moreover, since the 2-D Ricci scalar 
for the parametrization (13.11) is 



V^R = V^iRig*) - 20*p) (3.4) 

where □* refers to the D'Alembertian of the metric g^^, = e~'^''g^i,, it follows that after 
an integration by parts, the action contains terms at most linear in the time derivative of 
p. This in turn implies that the gravitational part of the action can formally be written 
in the form: 

Sg = j dtdx(^B^{<P,^,p,a,N)p + Bo{^,^,p,a,N)'^ . (3.5) 

Thus 

H, = 5i(0,0,p,a,iV) (3.6) 
which in principle at least can be solved for 0: 

= 0(H„p,0,a,iV) (3.7) 

Note that in the above and what follows we do not specify separately dependence on 
spatial derivatives of a field separately from dependence on the field itself. 
The Legendre transformation then yields a Hamiltonian of the form: 



Hg = J dx{U^<j){Up, p, 0, a, N) - 5o(0(H„ p, 0, a, iV), p, 0, a, N)). (3.8) 

where as indicated is an implicit function of the other phase space variables as obtained 
from (13. 6 p but does not depend on H,^. The crucial property of the above is that the p 
terms have cancelled and the only dependence on H^ is in the first term. 

We now make the above exphcit for the action Sg = Seh + Sgb- In particular, we 
have: 

Sg = Q j dxdt{Ai^^ + ^303 + ^20^ + Ai0 + Aq) (3.9) 
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Here 7 is defined as 7 = P i6(n-i)(Ti.-2) ^j^^ ^^^^ primes represent differentiation 
witli respect to tfie t and x coordinates respectively. 

For completeness we note that the conjugate momentum of is: 

= Ai + 2A2(j) + 3^30' + 4 A4(/)^ (3.15) 



3.2 Partial Gauge Fixing 

The above general discussion suggests a particularly convenient family of gauges, namely 
those for which the areal radius is time independent. We therefore consider the class of 
gauge fixing conditions: 

X/ = 0-/(x)^O (3.16) 



3.2 Partial Gauge Fixing 



If f{x) = x^/", then this gives x = r, i.e. the spatial coordinate is equal to the areal 
radius, but for the moment we will consider the more general case. 

The ~ denotes the fact that the gauge fixing condition can only be imposed on initial 
data. If it is a good gauge fixing condition then the requirement that it be preserved by 
the time evolution will yield a condition on one of the lagrange multipliers a or N. Xf 
must be second class with some linear combination of the constraints. Clearly: 
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From Eq.f l3.6p we can deduce that the required condition on the phase space variables is: 
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which determines the shift N in terms of a, p, lip and cj). The latter is no longer dynamical. 

In order to set x strongly equal to zero in the Hamiltonian (as opposed to the equations 
of motion) we must use the formalism of Dirac brackets. Fortunately, since is no longer 
a phase space variable, the Poisson bracket of all phase space variables with x is zero. 
This means that the Dirac bracket is the same as the Poisson bracket. 

The next step in the gauge fixing procedure is to use the diffeomorphism constraint 
to solve for II,;!, in terms of the other phase space variables. While this is indeed possible, 
it is not necessary in this case because the general form of the Hamiltonian (13. 8 p implies 
that when = 0, II,^ is not present in the Hamiltonian. Thus in this family of gauges, 
the gravitational part of the Hamiltonian is simply: 



Hg = - I dxBn. 
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where it is understood that in Bq, the shift is replaced by its value as determined by 

(EUD. 

To find the Hamiltonian with = we use fl3.19p . noting that -Bo(0 = 0) = Aq{p = 0) . 
This gives 
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where Hp is calculated to be 
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which is a cubic equation that can be solved for N/a. As noted in [TU] there can in 
principle be multiple real roots, depending on the parameters involved, which would 
yield different dynamical spacetimes. In that work, the boundary conditions effectively 
resolved this potential ambiguity. In the appendix we show that in our case as well there 
is a unique real branch, and hence unique dynamics, for the fiat slice, generalized P-G 
coordinates that we consider. 

We therefore assume that there is a unique solution that specifies implicitly N/ cr as a 
function of the remaining phase space variables and we use this to write the Hamiltonian 
in terms of one Lagrange multiplier, cr, and one constraint 
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where V is the coefficient of R in the action, equation (12.111) iV,^ = 7^0^ + 8(^rxy' 
and g^^ is given by 



After some algebra, the Hamiltonian can be written as 
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where we have defined the mass function Ai [1T| as: 
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Using: 
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one can verify that M. agrees with the mass function defined in equation (2.7) of [14J. 
It is therefore the EGB generalization of the Misner-Sharp mass function. It is invariant 
under coordinate transformations that preserve spherical symmetry and asymptotes to 
the ADM mass for asymptotically fiat spacetimes. Clearly Ai is constant for vacuum 
spacetimes. 



3.3 Adding the Matter Contribution 



Eq.f l3.25p is a quadratic that can be inverted to express (g^^) and (via fl3.23p ) N/a in 
terms of the mass function: 



l±Jl + 167(72 - l)20-2-2/n>f /nQ (3.27) 



From the above, we see that in vacuum, when = M is constant, the above yields the 
two branches of the well-known vacuum solution of EGB theory. For the following, we 
choose the plus-branch in f l3.27p in order that our solutions reduce to those in Einstein 
gravity in the limit 7 — )■ 0. However, nothing in principle precludes the other choice, 
which would have interesting consequences for the end-point of gravitational collapse. 
This in turn puts the ADM metric in the form: 

ds^ = e^" [-a^e^'{g^^)de ± 2a^l - e^P{g^^)dtdr + cir^) . (3.28) 

with {g^^) given in fl3.27p above. Thus, the partially gauge fixed geometry is now implic- 
itly specified in terms of p, lip and the the lapse a. 

Finally, noting that the condition for the existence of an apparent horizon is that the 
null expansion vanishes: 

|Vr|2 = (7ii(r')' = (3.29) 

it is easy to verify that by choosing the -f- (-) in front of the off-diagonal term fl3.28p one 
can guarantee that the slicings are regular across future (past) horizons. The exception to 
this is the gauge choice = (Schwarzschild gauge) which breaks down at all horizons. 

3.3 Adding the Matter Contribution 

In order to add the scalar field contribution to the action we dimensionally reduce equa- 
tion ([MD to get 

Sm = TT- / dxdt^/^(j)^\V^\^ = (3.30) 
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From this we can write the matter action as 
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From equation (13.321) it is clear that the total partially gauge-fixed Hamiltonian is 
given by 
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where (^) is understood to the be real branch of the solution to (13.211) . The Hamiltonian 
(I3.33P is the generalization to G-B gravity of the partially reduced Hamiltonian derived 
in [15]. Note that it is valid for any gauge fixing condition Xf the form (I3.16p . For 
concreteness we henceforth specify: 
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which chooses the areal radius as the spatial coordinate: x = r = </)^/". 

At this stage, the dynamical system consists of four phase space variables (p, Hp, xjj, H^) 
and there is a single first class constraint, namely the Hamiltonian constraint. We now 
look at the complete reduction to only the physical degrees of freedom associated with 
the scalar field. 



4 Completely Reduced Theory 

For our second gauge choice we follow [8] and choose 
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It is clear from (13. ip and (12. 5 p that this choices yields spatial slices that are fiat. It 
therefore generalizes P-G coordinates. To further understand the significance of this 
choice, note that using (I3.27P in (I3.28P and setting x to zero in the limit as 7 — > 
(choosing the plus sign), we get g^^ — j- 1 — 2GM./r. Inserting this into (I3.28P and setting 
^ to zero we find that the metric becomes 
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This shows that the gauge choice ^ = puts the metric into the Gauss-Bonnet, non- 
static generalization of the Painleve Gullstrand co-ordinates. Using ^ = also allows us 
to write the consistency condition for the first gauge choice, equation (I3.27p . as 
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5 Conclusion 



To ensure that the second gauge condition is conserved in time we must insist that 
{d/dt){e^p - 1) = {{e^P -1),H} = 0^ iS/6Ilp)H = 0. From flX^ we know that the Up 
dependence in the Hamiltonian is in the {N/a) terms. So we can write 
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where it is understood that we use (14.31) to write {N/a) in terms of the mass function 
M. 

Using X = ^ = we can now write down the fully gauge fixed Hamiltonian as 
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Using Hamilton's equations we find 
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These equations, along with the consistency conditions (14. 3 p and (14. 5 p and the Hamilto- 
nian constraint 
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determine the evolution of a collapsing scalar field. 



H^^' = 



(4.9) 



5 Conclusion 

Lovelock gravity has become increasingly popularized in recent years, as other theories 
have required higher dimensional analogues to General Relativity. We have presented 
here the Hamiltonian formulation of the lowest order non-trivial Lovelock theory, Gauss- 
Bonnet gravity, minimally coupled to a scalar field in the spherically symmetric case. 

The reduced Hamiltonian system (14. 51 - 14. 9]) is well suited for the study of gravitational 
collapse of a scalar field along with the resultant critical phenomena. Our calculations 
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utilize a gauge choice that selects non-static Painleve-GuUstrand coordinates, a coordi- 
nate system that is regular across the horizon during black hole formation. Moreover, this 
system is easily generalized to any matter field whose action does not contain derivatives 
of the metric. 

Future work will involve numerically studying gravitational collapse of a scalar field 
in this scenario to see how non-linearity in curvature modifies the well-known critical 
phenomena [1] in Einstein gravity. 
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7 Appendix 

We now solve equation 13.211 for N/a. Equation 13.211 can be written in the form 

c,(^)\c,[^)-n,^o (7.i; 

where 



Providing C1/C3 > 0, there is only one real solution, namely: 
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The condition that C1/C3 > can be simplified to 

^ 7n2</)2("-2)/" + ^\ - 1 > (7.4) 



470'^ V ^ - 1 

for the solution to equation 17.11 to be real. This won't be identically true so to proceed 
we must be careful with our gauge choices or possibly our allowed range of 7. For the 
fiat slice coordinates of Section (jl]) x = r (ie 0' = (r;,/2)0^~^/") and e^^ = 1 so that the 
condition is in fact identically satisfied. 
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